Superlattice Structures of Graphene based Nanoribbons 
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Based on first-principles calculations we predict that periodically repeated junctions of armchair 
graphene nanoribbons of different widths form superlattice structures. In these superlattice het- 
erostructures the width and the energy gap are modulated in real space and specific states are 
confined in certain segments. Orientation of constituent nanoribbons, their width and length, the 
symmetry of the junction are the structural parameters to engineer electronic properties of these 
quantum structures. Not only the size modulation, but also composition modulation, such as period- 
ically repeated, commensurate heterojunctions of BN and graphene honeycomb nanoribbons result 
in a multiple quantum well structure. We showed that these graphene based quantum structures 
can introduce novel concepts to design nanodevices. 
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Recent advances in materials growth and control tech- 
niques have made the synthesis of the isolated, two- 
dimensional (2D) honeycomb crystal of graphene possi- 
ble [3, 0, Q . Owing to its unusual electronic energy band 
structure leading the charge carriers resembling massless 
Dirac Fermions, graphene introduced novel concepts and 
initiated active research [J, H, H, 0, HI- F° r example, 
quasi ID graphene nanoribbons exhibit interesting size 
and geometry dependent electronic and magnetic prop- 
erties @, E3, EH [l2 - 

In this letter, based on extensive first-principles, as 
well as empirical tight-binding (ETB) calculations we 
showed that periodically repeating heterojunctions made 
of graphene nanoribbons of different widths can form 
stable superlattice structures. These superlattices are 
unique, since both size and the energy band gap are mod- 
ulated in real space. As a consequence, in addition to the 
propagating states, specific states are confined in certain 
regions. Confinement increases and turns to a complete 
Mott's localization when the extent of the nanoribbons 
with different widths increases. Widths, lengths, chiral- 
ity of constituent nanoribbons and the symmetry of the 
junction provide variables to engineer quantum struc- 
tures with novel electronic and transport properties. We 
also showed that multiple quantum well structures can 
be formed through composition modulation of commen- 
surate heterojunction of BN honeycomb and graphene 
nanoribbons both having the same width. In particu- 
lar, our transport calculations demonstrate that specific 
finite size quantum structures can operate as resonant 
tunnelling double barrier (RTDB) with resonances, char- 
acteristics of quantum well states. Since nanoribbons 
having width less than ~10 nm with uniform and pat- 
terned edges have been produced [l3[ , present results re- 
lated with graphene based nanoribbons and their hetero- 
junctions of diverse geometry can, in fact, be a candidate 
for a class of nanodevices with a richness of novel prop- 
erties. 

We have performed first-principles plane wave calcu- 
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lations [Tj, [l5| within density functional theory (DFT) 
0, E3 using projector augmented wave (PAW) po- 
tentials [III]. The exchange correlation potential has 
been approximated by Generalized Gradient Approxi- 
mation (GGA) using PW91 functional[19j both for spin- 
polarized and spin-unpolarized cases. All structures have 
been treated within supercell geometry using the peri- 
odic boundary conditions. A plane-wave basis set with 
kinetic energy cutoff of 500 eV has been used. In the 
self-consistent potential and total energy calculations the 
Brillouin zone (BZ) is sampled by (1x1x35) special It- 
points for nanoribbons. This sampling is scaled according 
to the size of superlattices. All atomic positions and lat- 
tice constants are optimized by using the conjugate gra- 
dient method where total energy and atomic forces are 
minimized. The convergence for energy is chosen as 10 -5 
eV between two steps, and the maximum force allowed 
on each atom is 0.02 eV/A. 

Here we consider the bare and hydrogen-terminated 
armchair graphene nanoribbons, AGNR(n), n being the 
number of carbon atoms in the primitive unit cell. All 
AGNR(n) are semiconductors and their band gap, E g , 
vary with n. It is small for 3n/2 — 1, but from 3n/2 to 
3n/2 + 1 it increases and passing through a maximum it 
becomes again small at the next minimum correspond- 
ing to 3n/2 + 2 [T3|. As E g oscillates with n its value 
decreases eventually to zero as n — > oo. AGNR(n) is non- 
magnetic. In Fig. [T] (a) and (b) the band gap values of 
bare AGNR(IO) and AGNR(14) are 0.44 and 1.10 eV, re- 
spectively. According to isosurface plots one distinguish 
the relatively uniform states (having their charge uni- 
formly distributed across the nanoribbon) from the edge 
states (having their charge accumulated at both edges of 
the ribbon) . Owing to the interaction between two edges 
of the narrow ribbon, the bands of edge states split. Upon 
termination of the carbon dangling bonds by hydrogen 
atom, the edge states of AGNR(IO) and AGNR(14) dis- 
appear, and their energy band gaps change to 0.39 eV 
and 1.57 eV, respectively. Since hydrogen terminated 
AGNR's are stable at ambient conditions, in the rest of 
the paper we consider their superlattices unless stated 
otherwise. 
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FIG. 1: (a) Bare and hydrogen terminated AGNR(IO): 
Atomic geometry, electronic band structure and isosurface 
charge densities of edge and "uniform" states. The primitive 
unit cell is delineated with dashed lines and includes n=10 
carbon atoms. Carbon and hydrogen atoms are shown by 
large and small balls, (b) Same for AGNR(14). (c) Variation 
of band gaps E g of AGNR, with the tensile strain, e. (d) 
Variation of the total energy of AGNR with e and its second 
derivative with respect to e, i.e. k . All data in this figure are 
calculated by first-principles method. 



The rate of change of the band gap with strain, e in 
Fig. [He), i.e. dEg/de is significant and its sign changes 
with n. The mechanical properties of AGNR(n) are also 
important for any device application. Here we calculated 
the variation of total energy with strain, e and its second 
derivative with respect to strain, k — 8 2 Et /de 2 . Results 
summarized in Fig.[T](d) indicate that AGNR's are quasi 
ID, stiff materials. For example, n =697 eV/cell for 
AGNR(14). These values can be compared with n =127 
eV/cell calculated for linear carbon chain [28| and con- 
firm the robustness of heterostructures. 

The remarkable properties of graphene ribbons dis- 
cussed above, especially their band gap rapidly vary- 
ing with their width, suggest that the heterostructures 
formed by the segments of AGNR's with different n may 
have interesting functions. Here the crucial issues to be 
addressed are how the electronic structure will be affected 
from the discontinuity of crystal potential at the junc- 



tion, what the character of the band discontinuity and 
the resulting band alignment will be, and whether these 
discontinuities will result in confined states. To this end 
we first consider a superlattice AGSL(ni,n2;si,S2) made 
by the segments of AGNR(ni) and AGNR(n 2 ). Here, 
Si and S2 specify the length of segments in terms of the 
numbers of their unit cell. Fig. [2] shows the superlattice 
AGSL(10,14;3,3). It is made by periodically deleting one 
row of carbon atoms at both edges of AGNR(14) to form 
periodic AGNR(10)/AGNR(14) junction. Upon junction 
formation, dramatic changes occur in the band structure 
of this superlattice. While highest (lowest) valance (con- 
duction) band is dispersive and their states propagate 
across the superlattice, the second valence and conduc- 
tion bands are flat and their states are confined to the 
wider part of AGSL(10,14;3,3) consisting of a segment of 
AGNR(14). These latter flat band states are identified 
as confined states. 

The superlattice AGSL(10,14;3,3) in Fig. E^a) can be 
viewed as if a thin slab with periodically modulated width 
in the sy-plane. The electronic potential in this slab is 
lower (V <0) than outside vacuum (V = 0) and thus is 
modulated in real space. Normally, electrons in this thin 
potential slab propagate along the x— axis. If the junc- 
tion between wide segment [AGNR(14)] and narrow seg- 
ment [AGNR(IO)] is sharp, some electronic states in the 
wide segment happen to reflect from sharp discontinuities 
at the junction and eventually become confined. The itin- 
erant states are one dimensional, and the band gap and 
width of this superlattice are modulated in real space. 
On the other hand, the atomic arrangement and lattice 
constants at both sides of the junction are identical; the 
hetero character concerns only the different widths of the 
ribbons at different sides. According to this simple pic- 
ture, the larger An = n 2 — ni (and also the larger si and 
S2, or the sharper the junction), the stronger becomes the 
confinement. As justified in further discussions, these 
arguments are relevant for superlattices with long con- 
stituent segments having lengths larger than de Broglie 
wavelengths. Note that AGNR's can be constructed by 
using two different unit cells (i.e. those consisting of com- 
plete hexagons or incomplete hexagons). For the wider 
parts of the superlattice the one with complete hexagons 
is preferred in order to prevent a lonely carbon atom 
at the interface. For AGSL superlattices with reflection 
symmetry the narrow region is made by unit cells having 
complete or incomplete hexagons depending on whether 
(n2 — ni)/ A is an even or odd number, respectively. 

In addition to n\, ri2, si, S2, the symmetry of the 
junction, An = ri2 — ni, even-odd disparity of ni/2 and 
712/2, the type of the interface between two different rib- 
bons and the overall shape of the superlattice (namely 
sharp rectangular, smooth wavy, or sawtooth like) in- 
fluence the electronic properties. As shown in Fig. [2] 
(d), as S1—S2 increases, various superlattice band gaps 
decrease. In terms of the weight of the states in seg- 
ment S2, i.e. f \^>(r)\ 2 dr, we see that the confine- 
ment of states increase with increasing S1—S2', but it 
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FIG. 2: (a) Atomic structure of AGSL(m = 10, n 2 = 14; si = 
3, S2 = 3). The superlattice unit cell and primitive unit cell 
of each segment are delineated, (b) Band structure with flat 
bands corresponding to confined states, (c) Isosurface charge 
density of propagating and confined states, (d) Variation of 
various superlattice gaps with si = S2. (e) Confinement of 
states versus si = S2- 



disappears for si=S2=l. Differing of s% from s 2 causes 
Eg and the band structure to change. For example, 
AGSL(10,14;si,s 2 ) has E g =0.66, 0.48, 0.38, 0.32, 0.27 
eV for s 2 =3, but si=3, 4, 5, 6, 7, respectively. Con- 
versely, E g =0.72, 079, 0.83, 0.84 eV for si=3, but s 2 =4, 
5, 6, 7, respectively. On the other hand, the energy of 
the flat-band states confined to s 2 and their weight are 
practically independent of s±. Variation of An = ni — n 2 
results in a wide variety of electronic structures. For 
example, in contrast to AGSL(10,14;3,3) the highest va- 
lence and lowest conduction bands of AGSL(10,18;3,3) 
are flat bands with £^=0.70 eV; dispersive bands occur 
as second valence and conduction band having a gap of 
1.19 eV between them. The situation is even more com- 
plex for AGSL(10,22;3,3): While the first and the sec- 
ond bands are dispersive in both valence and conduction 
bands with a direct band gap 0.23 eV, flat bands occur 
as third bands with a gap of 2.61 eV between them. The 
states associated with the flat bands near the Fermi level 
are confined at the wider part of the superlattice. It 
should be noted that the band gap of nanoribbons, E g 



is underestimated by the DFT calculations 20] . Since we 
consider structures which have already a band gap, its 
actual value does not affect our discussion in any essen- 
tial manner, but an enhancement of predicted properties 
(such as the extent of confinement) can be expected. 

The above trends corresponding to small n\ and n 2 
become even more interesting when ri\ and n 2 increase. 
The electronic band structure of AGSL(ni,n 2 ;si,s 2 ) with 
ni=42 or 82, but n 2 > n\ and si=s 2 > 3 calculated 
using ETB method shows that for small An, confine- 
ment is weak and bands are dispersive, but confine- 
ment increases as An increases. Interestingly, E g of 
AGSL(ni=42,n 2 ;3,3) is, respectively, 0.46, 0.12, 0.49 and 
0.04 eV for n 2 =46, 50, 54 and 58. In Fig. [3] we demon- 
strate that the confined states can occur not only in nar- 
row (small ni and n 2 ), but also in wide superlattices 
having significant modulation of the width. In ETB 
method used here, the Bloch states having band index 
I and wavevector k are expressed in terms of the linear 
combination of the orthonormalized Bloch sums Xi(k, r) 
constructed for each atomic orbital p z localized at differ- 
ent carbon atoms, i, with the proper phase of k, namely 
V>z( k > r ) = ^^^(kjx^k.r). Accordingly, the contribu- 
tion of the orbital at site i to the normalized charge den- 
sity of Vi(k,r) is given by p t = J2i K*( k )| 2 - In Fig. G2 
Pi is scaled with the radius of circles located at atomic 
site i. 

Electronic structure is also strongly dependent on 
whether the geometry of superlattice is symmetric (hav- 
ing a reflection symmetry with respect to the superlat- 
tice axis along the x-direction) , or saddle (one side is 
straight, other side is periodically caved) all having the 
same An. While the saddle structure of AGSL(10,18;3,3) 
has largest direct gap between dispersive conduction and 
valence bands, its symmetric structure has smallest gap, 
but largest number of confined states. Horn- like smooth 
connection between wide and narrow segments may give 
rise to adiabatic electron transport and focused electron 
emission 
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The confined states have been treated earlier in com- 
mensurate or pseudomorphic junctions of two different 
semiconductors, which form a periodically repeating 2D 
superlattice structure [221] . These superlattices have 
grown layer by layer and formed a sharp lattice matched 
interface. Owing to the band discontinuities at the inter- 
face, they behave as if a multiple quantum well structure 
according to the Effective Mass Theory [23[ . Two dimen- 
sional conduction band electrons (valence band holes) 
confined to the wells have displayed a number of elec- 
tronic and optical properties [24], l25l l26l . [27j ■ In Fig. [U 
we present a ID analog of the 2D semiconductor su- 
perlattices through compositional modulated nanorib- 
bons. BN honeycomb ribbon and graphene ribbon of 
the same width are lattice matched and can form super- 
lattices with multiple quantum well structure with con- 
fined states. We considered a periodic junction of the 
segment of armchair BN ribbon with ni = 18 and si=3 to 
the segment of armchair graphene with n 2 =18 and ,s 2 =3 
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FIG. 3: Energy band structure of the AGSL(34,70;11,9) su- 
perlattice and the charge densities of selected bands. As seen 
clearly, states associated with flat bands 1,3, and 4 are con- 
fined but the state with dispersive band indicated by 2 is 
propagating. Calculations have been performed using ETB 
method. 



to form a superlattice structure. While periodic BN and 
graphene ribbons by themselves have band gaps of ~5 eV 
and 0.8 eV, respectively, the band gap of BN/AGNR(18) 
is only 0.8 eV indicating a normal (type-I) band align- 
ment. Under these circumstances, a state propagating 
in one segment becomes confined if its energy coincides 
with the band gap of the adjacent BN segment. In 
Fig. [4] (b) , the dispersive minibands and non-dispersive 
quantum well states are clearly seen. That these quan- 
tum well states are confined in the graphene zone (which 
has small band gap as compared to the band gap of BN 
ribbon) are demonstrated by isosurface plots of charge 
densities. In contrast to the confined states, the propa- 
gating states have charge densities in both graphene and 
BN zones of the superlattice. This is another class of 
heterostructure obtained from graphene based ribbons 
and their functions can even be advanced by implement- 
ing the size modulation in addition to the compositional 
one. 

Finally, we present an interesting device that can op- 
erate as a RTDB made preferably from a segment of 
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FIG. 4: One-dimensional superlattice structure formed from 
the junction of BN and graphene armchair nanoribbons. (a) 
Atomic structure and superlattice parameters, (b) Band 
structures of constituent BN and graphene armchair nanorib- 
bons having 18 atoms in their unit cells and the band struc- 
ture of the superlattice BN(18)/AGNR(18) each segment hav- 
ing 3 unit cells (si = S2 =3). (c) Energy band diagram in real 
space forming multiple quantum wells, QW, in graphene seg- 
ments (zones). Isosurface charge densities of states confined 
to QW's and propagating states are presented for selected 
bands. 



AGSL(10,18;4,4). Such a device is relevant for appli- 
cations and uses LUMO and HOMO states confined in 
the wide region. For the sake of illustration we con- 
sidered fictitious metallic electrodes of two widely sepa- 
rated (weakly coupled) carbon chains 28| altogether hav- 
ing 4 quantum conductance channels. Six principal lay- 
ers of electrodes are included at both sides of RTDB 
as parts of the central device. Metallic electrodes make 
perfect contacts with the central RTDB device. Trans- 
mission coefficient, T calculated self-consistently using 
non-equilibrium Green's Function method reflects the 
combined electronic structure of central RTDB device, 
electrodes and their contacts. Double-^ plus polarization 
numerical orbitals have been used as the basis set and 
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atomic structures are further optimized before transport 
calculations. Fig. [5] presents the calculated transmission 
curve. Confined LUMO and HOMO states and other 
confined states identified through the energy level di- 
agram and isosurface charge density plots give rise to 
sharp peaks originating from resonant tunnelling effect. 
States extending to whole RTDB are coupled with the 
states of electrode are shifted and contributed broad 
structures. 



In conclusion, various types of quantum structures 
have been revealed through size and composition modu- 
lation of graphene based nanoribbons. The confinement 
of electron and hole states are shown to lead to interest- 
ing effects for new device applications based on graphene. 
Several structural and compositional parameters to mon- 
itor electronic and transport properties are expected to 
make these quantum structures an attractive field of re- 
search. 



FIG. 5: (a) Atomic structure of a resonant tunnelling 
double barrier (RTDB) device made from a segment of 
AGSL(10,18;4,4) and connected to two metallic electrodes, 
(b) Transmission coefficient, T versus energy calculated un- 
der zero bias, (c) Energy level diagram of the RTDB and 
charge densities of confined and extended states. 
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